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Abstract 

We derive spectral sum rules for a system with two quarks coupled to an imaginary 
isospin chemical potential in the e regime. The sum rules show an explicit dependence 
on the pion decay constant which should make it possible to measure F n from the 
eigenvalue spectrum of this particular Dirac operator. 
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1. The determination of low energy constants of QCD such as the pion decay con- 
stant F n remains an important problem. In particular the computation of such 
quantities from lattice calculations is notoriously difficult due to the exceeding 
computational challenge posed by simulations of small quark masses. In addition 
when approaching the chiral limit on the lattice, finite size effects inevitably be- 
come more and more significant. Hence, a method which takes finite size scaling 
explicitly into account seems to be highly profitable. Such a technique is provided 
by the so called e-regime of QCD [1]. This regime applies in a region where the 
Compton wave length of the pion l/m n is larger than the one dimensional size 
L of the physical volume V = L 4 , while still being much smaller than the typi- 
cal hadronic scale Aqcd-, i-e. l/jro^ > L ^> 1/Aqcd- The lowest order effective 
partition function of the e-regime is known analytically. The fact that it depends 
explicitly on the infinite volume chiral condensate E = {ipip) can be exploited to 
determine this constant with high accuracy: back in 1992 Leutwyler and Smilga 
derived a set of spectral sum rules for the Dirac operator by restricting the parti- 
tion function to sectors with fixed topological charge v [2]. In this way the chiral 
condensate is linked to the spectrum of the Dirac operator in finite volume which 
can be determined from lattice simulations [3,4]. 

However, with the standard Dirac operator other low energy constants such 
as the pion decay constant F n appear only in higher order corrections [5]. As 
a consequence, the computation of F n was believed to be much more demanding 
[4,6,7]. It is therefore still more widespread to extrapolate the pion decay constant 
down to the chiral limit from simulations in the p regime of chiral perturbation 
theory [8,9]. 

Recently a new approach has been proposed [10-12] which avoids these dif- 
ficulties: if the quarks are coupled to an external source, dependency on F n ap- 
pears already at the lowest order in the partition function [5,13]. The authors of 
Refs. [10,11] have used this fact to derive correlation functions of the eigenvalue 
densities which are sensitive to F n for both quenched and unquenched chiral per- 
turbation theory. Here, we will use the same partition function to deduce a set 
of sum rules in a way analogous to Ref. [2]. These rules depend then likewise 
on F n and make it possible to determine this important quantity from a lattice 
simulations in finite volume. 

2. We consider a system with two quark flavours u and d which are coupled to an 
external source Hi SQ . The source can be interpreted as an imaginary isospin chem- 
ical potential which couples differently to the two flavours or as twisted boundary 
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conditions [12,14]. This gives rise to two independent eigenvalue equations 

(D+ + m u ) ip +n = (lp[A) + ifiiso^o + rn u ) ip +n = (i\ +n + m u ) (1) 
(D- + m d ) = (#[A) - ifii SO -f + m d ) ip- n = («A_„ + m d ) ?/>_„, (2) 

where A denotes the gauge field. The advantage of an imaginary isospin chemical 
potential is twofold. Firstly, an isospin chemical potential preserves the positivity 
of the the fermion matrix [14]. The system can thus be simulated on the lat- 
tice with the usual Monte Carlo techniques without running into sign problems. 
Secondly, if it is chosen to be imaginary the massless operators D + and D_ are 
anti-hermitian and the eigenvalues A± lie on the real axis. 

In the case of degenerate masses m u = m d the two flavours can be converted 
into each other by the transformation \ii so to —fii, so . Here we use a non mass- 
degenerate quark pair for purely technical reasons: as we will see below this gives 
us the only handle to distinguish between the two sets of eigenvalues. Thus keeping 
the masses distinct allows us to extract sum rules for only one set of eigenvalues, 
A + or A_ accordingly. 

For both sets of eigenvalues the non zero modes come in positive and negative 
pairs, ±A +n and ±A_ n respectively. Hence in a sector of fixed topology, the 
partition function can be written as 

Z u (m u ,m d ) = (m v u m» d H(\ 2 +n + m 2 u ) l[(X 2 _ n + m 2 d )\ , (3) 

\ n n I „ 

where {■)„ denotes the gauge average over all configurations with topological 
charge v and the products are restricted to strictly positive values of A. 

Our derivation of the sum rules deviates from the procedure originally used 
by Leutwyler and Smilga in Ref . [2] , we will follow the somewhat easier approach 
presented in Ref. [15] instead. Massive spectral sum rules for A + (or A_) can be 
derived from formula (jSJ) by taking logarithmic derivatives with respect to the 
masses m u (or m d ) . A first order sum rule for instance is given by 

As in the products of Eq. (j3J) the sum on the rhs runs over the positive eigenvalues 
only. 

On the other hand in the e expansion of chiral perturbation theory the par- 
tition function can be shown [11, 13] to be 

Z v (x v , x d ) = e~ 2VF ^o f ds e 2VF^ s y s Iv ( SXd ) Iv ( sx ^ ( 5 ) 

Jo 
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where we introduced the scaling variables of the e-regime Xi = VY>rrii, i — u,d and 
I u are modified Bessel functions. The dependence on F n is through the product 
VF 2 fj, 2 S0 only. In particular a sign change in \i iso leaves the partition function 
invariant. The two flavours are indeed only distinguishable through their masses. 
Upon inserting this explicit formula for Z v into Eq. (jljl. we obtain the first order 
sum rule 

1 \ V 2 Y? ft ds e 2VF ^s° s2 s 2 I v (sx d ) Wfxj 
r A +«+</„ 2x « j; ds e 2VF ^ 2 s I u (sx d ) I u (sx u ) ' 

This formula simplifies a good deal if we take both quark masses to zero. The 
Bessel functions disappear up to remnant powers of the parameter s 



V 2 Y? fUse 2VF ^ 2 s ^ 

( 7 ) 



S \* I 4(„ + i) ri ds e 2^>? soS2s2 , + i 

For general z/ this expression can only be evaluated in terms of incomplete and 
ordinary T-functions 



\2 



VT 2 T(v + 2) -T{v + 2,-2VF 2 fi 



2, ,2 N 

iso) 



\%J 8(1/ + r(i/ + i) - r(i/ + i, -2vjm 



2 ) 

ISO/ 



(8) 



However, for a given topological charge it is straightforward to calculate the pa- 
rameter integrals. The rule reduces then to a rather simple expression, where F n 
appears only in polynomials and exponentials and which can easily be fitted to 
lattice data. Let us illustrate this for the case of vanishing topological charge, 
there we simply have 



\2 



US 2 l + e 2VF ^(2VF 2 tf so -l) (q) 

Eqs. (J2J) and (jHJ) are the direct equivalents of the first order sum rule given in 
Ref . [2] , but evaluated for a system with two quarks coupled the chemical potential 
\i iso . Indeed, in the limit where fi iso vanishes Eq. ((7|) becomes 



\2 



V 2 Y? tidss 2u+3 V 2 T? 



Xi n / 4(v + l)fidss2"+i 4(i/ + 2) 
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which is precisely the Leutwyler-Smilga result for Nf = 2 flavours. 

Due to the symmetry of the partition function under an exchange of the quark 
masses, the corresponding sum rule for A_ can be obtained from Eq. © by simply 
substituting x u for If we had treated the quarks as mass degenerate from the 
beginning we would have arrived at the sum of those two sum rules, which is just 
twice Eq. (JOJ). 

Higher order sum rules 

We can derive two different types of second order sum rules by either taking the 
second derivative with respect to m u , In Z v , or a mixed derivative dm dm In Z u . 
Let us start with the second possibility. Applied to Eq. (j3J) the mixed derivative 
yields the subtracted correlation of the two sets of eigenvalues 




Since we have calculated the disconnected parts already, the only new contribution 
from Eq. (jllj) is a sum rule for 




In the limit of vanishing quark masses we can again expand the Bessel functions 
and obtain a much simpler expression for the massless mixed second order sum 
rule 

// V _L_\ / V _L_\\ ^dse^^V^ 5 _ 

\ V n A +J V » A " J /„ _ 16 ^ + ^ 2 Jo ds e 2VF ^°° s2 S 2 » +1 ~ 

V 2 ^ T(u + 3) - I> + 3, -2VF^j ) 
64F>f s > + l) 2 T(u + 1) - T(y + 1, -2VF^L) ' 

(13) 

Apart from the coefficients, the difference to the massless first order rule is given 
by the higher power of the parameter s in the first line of Eq. (|13p. These powers 
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are indeed a distinctive feature for any massless sum rule of a given order (third 
order sum rules for instance carry a power of s 2u+7 ). Concerning the evaluation 
of the integral we can make exactly the same remarks as for the first order rule, 
they are easily calculated in a fixed topological sector. For completeness we give 
again the result at v — 



V 2 ^ (1 - 2VF 2 p 2 so (l + VF^j ))e^lo 1 
XlnJ \^ A 2 „ j I n 32FXo e 2 ^>L - 1 

(14) 

Again, in the limit /^j so — ► where A + and A_ become degenerate the formula 
reproduces the result of Ref. [2] In addition, this second order sum rule can be 
compared to the results of Ref. [11]. There the mixed two point spectral correlation 
function 



p {2) {Xi,X2,m u ,m d ,ifi iso ) = ( ^SjXi - A +n )^5(A 2 - A_ 

\ n I 



(15) 
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is derived with the replica method. By integrating out Ai and A2 in 

j\ f a\ P (2) (X l ,X 2 ,m U} m d ,ip lso ) 
dXi J dX 2 ^2 ( 16 ) 

we should reproduce our subtracted mixed sum rule of Eq. (JTTJ). Indeed, we do 
find exact agreement in the massless limit. 

With the first choice above In Z v , we have means to extract the correlation 
between eigenvalues of one series A + (or equivalently A_). More precisely, from 
taking the second derivative with respect to m u and subtracting the zero mode 
contributions and the known first order terms, we can deduce a formula for 




X\ n + ml) ^ \X\ n + ml 

_ V 4 S 4 /„ ds e 2VF ^Q s2 s 2 I v {sx d ) [ s x u I v {sx u ) - 2{v + l)I u+1 (sx u )} 
~ 4 *l Jo ds e 2VF ^io* 2 s I u (sx d )I u (sx u ) 

Again the formula simplifies drastically in the massless limit where we get 

V 4 E 4 /„ ds e 2VF ^ 2 s 2,+5 
= 16(i/ +l)(v + 2) Jl ds e 2VF t&o° 2 s 2.+i 
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(17) 




Let us remark here that at finite /ii so we cannot extract a separate sum rule for 
the fourth order term (^2 l/A^)^ from the partition function in Eq. (0). Since the 
u and the d quark obey two different eigenvalue equations, each of them behaves 
effectively like a one flavour system. In order to isolate the fourth order term 
above we would need another flavour with the same eigenvalues A +n but different 
mass rh u . At vanishing isospin chemical potential however, where A + and A_ 
become degenerate a fourth order term can be calculated from the difference of 
the two second order sum rules described here and it again coincides with the 
result of Ref. [2]. 

It is in principle possible to continue this series and to calculate sum rules 
of any order by taking higher and higher derivatives, and subtracting the known 
lower order terms. The computation however becomes tedious and its usefulness 
doubtful. 

3. In this paper we have derived a set of spectral sum rules for a system of 
two quarks coupled to an imaginary isospin chemical potential from the finite 
volume partition function. This can be seen as an extension of the pioneering 
work of Leutwyler and Smilga in Ref. [2] . The sum rules derived here inherit the 
F n dependence from the partition function and provide means to determine this 
important constant from the spectrum of the Dirac operator introduced here. 
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